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ABSTRACT 

We present the result of the observation of the ultrafaint dwarf galaxy Willman 1 performed with the 
17 m MAGIC telescope during 15.5 hr between March and May 2008. No significant 7-ray emission 
was found. We derived upper limits of the order of 10 -12 ph cm~ 2 s _1 on the integral flux above 100 
GeV, which we compare with predictions from several of the established neutralino benchmark models 
in the mSUGRA parameter space. The neutralino annihilation spectra are defined after including the 
recently quantified contribution of internal bremsstrahlung from the virtual sparticles that mediate the 
annihilation. Flux boost factors of three orders of magnitude are required even in the most optimistic 
scenario to match our upper limits. However, uncertainties in the dark matter intrinsic and extrinsic 
properties (e.g., presence of substructures, Sommerfeld effect) may significantly reduce this gap. 
Subject headings: MAGIC - Willman 1 - SDSSJ1049+5103 - Dwarf Spheroidal Galaxies 
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1. INTRODUCTION 

Dwarf spheroidal galaxies (dSphs) are believed to 
be the smallest (size 1 ~ kpc), faintest (luminosi- 
ties 10 2 — 10 8 L©) astronomical objects whose dynam- 
ics a re dominated by dark matter (DM; iGilmore et al.l 
l2008l and references therein). They are found as satel- 
lites orbiting in the gravitational field of a larger host 
galaxy (e.g., the Milky Way (MW)). Their member stars 
show large circular velocities and velocity dispersions 
that, combined with their modest spatial extent, can 
be interpreted with the presence of a large DM halo 
of the order of 10 5 — 10 8 M . In recent years, the 
Sloan Digital Sky Survey (SDSS; lYork et al.l 120001) led 
to the discovery of a new population of MW satel- 
lites, comprising about as many new objects as were 
previously known (iBelokurov et al. 2004; Will man et al. 
I2005al IZucker et alJl2006t llrwin et al . 2007t iWalsh et al.l 
2007). This population of extremely low-luminosity 
galaxies is very interesting for DM searches and to study 
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the galaxy formation at the lowest mass scale. The ex- 
istence of a new class of ultrafaint MW satellites is also 
relevant because it provides a part ial solution for the 
so-called missi ng sat e llite problem (iKlypin et al.l [l999t 
Simon fc Gehal 120071 : iStrigari et all 120071 : iMadau et all 



20081 ) by partially filling the gap between the predicted 
and the measured number of galactic subhalos. In 
this fam ily is Willman 1, discovered by Wil lman et al.l 
(|2005an and soon established as potenti ally the most 



DM-dominated dSph satellite of the MW (IMartin et al 
2007t ISanchez-Salcedo fc Hernandedl2007t IStrigari et al 



2008|). 

The physics of DM has gathered much interest in 
recent years, in p articular after Wilki n son Microwave 
Anis otropy Probe (Sp ergel et al.l l2007t iKomatsu et al.l 
2008) measured its relic density with great precision. It 
is generally believed that DM manifests itself as a general 
class of weakly interacting massive particles that includes 
several candidates which satisfy both exp erimental con- 
straint s and theoretical requirements (see iBertone et al.1 
l2005al and references therein). Among them, one of 
the best theoretically motivated, for whom the relic den- 
sity is calculated without fine tuning from its nature, is 
the neutralino, arising in SuperS ymmetric (SUSY) the- 
ories beyond the Stan dard Model (|Wess fc Zuminolll974l ; 
lHaber fcTk anc. 19 81) and in particul ar in the mSUGRA 
extension (jChamseddine et al.l 1 19821 ). The mSUGRA 
neutralino annihilations can be observed through the 
production of 7-rays. The main emission comes from 
secondary products of hadronization processes and from 
final state radiation. In addition, line emissions are 
found through direct processes such as XX ~^ 11 1 an d 
XX which provide 7-rays of energies E = m x and 

E = m x — m 2 z o I m x respectively, even if those processes 
are at the l oop level and ther e fore st rongly suppressed. 

Recently IBringmann et all (l200 8a) , following an ear- 
lier idea from iBergstroml "( 1989h . showed that in some 
regions of the mSUGRA parameter space, a hitherto 
neglected contribution to 7-ray emission comes directly 
from charged sparticles mediating the annihilation into 
leptons, in processes like XX ~* l + l~l- They defined this 
intermediate state radiation as internal bremsstrahlung 
(IB). The IB mechanism permits to restore the helicity 
balance in processes that would otherwise be strongly for- 
bidden, and supplies 7-rays toward high energies (E > 
0.6 m x ) of up to several orders of magnitude. This 
boosted emission is particularly interesting for ground- 
based 7-ray observation, as performed by Imaging Atmo- 
spheric Cherenkov Telescopes (IACTs), which are sensi- 
tive above ~ 100 GeV where normally the IB boost takes 
place. 

As the 7-ray flux is proportional to the square of 
the DM density, IACTs focus on concentrated DM ob- 
jects. Dwarf galaxy satellites of the MW represent very 
good candidates, since they are the most DM domi- 
nated systems known in the universe, with very high 
mass-to- light ratios (M/L), close distance and reduced 
7-ray background from unresolved conventional Galac- 
tic sources (i.e. stellar evolutionary remnants). Some 
dSphs h ave already been st udied in 7-ray s: Draco by 
MAG IC ([Albert et al 1 l2008al) and Whipple (IWood et all 
l2008h : UMi by Whipple (IWood est alJ l2008h and Sagit- 
tarius by HESS (|Aharonian et al.ll2007D without any sig- 



nificant observation of DM annihilations and only flux 
upper limits were estimated. 

In this paper, we report results of the observation of the 
sky region around Willman 1 performed by the MAGIC 
telescope for a total of 15.5 hr between March and May 
2008. After a brief description of Willman 1 in Section [2j 
in Section |31 we estimate the flux using benchmark mod- 
els for the neutralino and a typical DM density profile. In 
Section 21 we briefly describe the MAGIC telescope and 
the Willman 1 data sample. In SectionJSJ we present and 
discuss the results of the observation and we set upper 
limits for the flux. In Section [6j we report our conclu- 
sions. 

2. WILLMAN 1 

In 2004, Willman et al. discovered a new MW com- 
panion SPSS J1049+5103 (10' l 40 m 22.3 s ,51°03'03.6"; 
Will man et all l2005al) as a faint overdensity of red, 
resol ved stars, which was observed again the next 
year (| Willman et alll2~005bl ) and named Willman 1. At 
that moment, this object represented the 10th dSph of 
the MW and, the first one discovered in 10 years. Fur- 
ther observations performed with the K eck/DEIMOS 
telesc ope confirmed the SDSS results (jMartin et all 
|2007[). while a mor e recent observation is reported 
bv lSiegel et all rf2008lh 

Willman 1 is located at a distance of 38 ± 7 kpc in 
the Ursa Major constellation. It is characterized by a 
very low number of resolved stars, the total luminosity 
being L = 855 L Q , and a very small half-light radius of 
ri/2 = 21 ± 7 pc, almost 2 orders of magnitude smaller 
than other known dSphs. The source was defined by 
previous authors as an "extreme" dwarf galaxy, because 
some of its characteristics lie between those typical for a 
globular cluster (GC) and those expected in an extremely 
faint dSph. The large spread in metallicity of its stars 
favors the dSph interpretation rather than that of a GC, 
whi ch would contain s tars of a similar age and metallic- 
ity ()Martin et all2007l) , even if this evidence was recently 
put under discussion by ISiegel et al.l ()2008f ). Willman 1 
is the least massive satellite galaxy known to date, with 
a total mass (M ~ 5 x 1O 5 M ) about an order of mag- 
nitude smaller than those of the least massive satellite 
galaxies previously known. However, the corresponding 
mass-to-light ratio, M/L~ 500 — 700 Mq/Lq, is one of 
the highest of all dSphs. This makes Willman 1 one 
of the most attra ctive dSph galaxies to lo ok for DM at 
present (see, e.g. . IBringmann et alll2008~bl ). its predicted 
DM annihilation flux being at least a factor of 3 larger 
than t he second best DM candidate, according to recent 
work (|Strigari et alll2007f ). 

3. THEORETICAL MODELING OF THE GAMMA-RAY 
EMISSION FROM WILLMAN 1 

The 7-ray flux originating from DM particle annihi- 
lations can be factorized into a contribution called the 
astrophysical factor J( V P) related to the morphology of 
the emission region and a contribution called the particle 
physics factor $ pp depending on the candidate particle 
characteristics: 



$(>£ ) = J(*)-$ pp (>£ ) 



(1) 



where Eg is the energy threshold of the detector and 'J 
the angle under which the observation is performed. 
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3.1. Astrophysical Factor 

At present, a concise and exact characterization of the 
DM density profile of Willman 1 is a delicate task, since 
observational da ta are still scar c e. Th is paper is based on 
the studies from IStrigari et al.1 (|2008|) , who modeled the 
profile by using only 47 stars after removing those with 
unclear kinematics. Furthermore, to avoid membership 
problems, only the observational data related to the inner 
half of the galaxy were taken into account. It is impor- 
tant to note that a null/insignificant tidal stripping was 
assumed in order to carry out the modeling of the DM 
distribution as a system in dynamical equ ilibrium, a fact 
which is still under debate. For example, I Willman et al.1 

(2006) claim the existence of strong tidal debris, the evo- 
lution of the dwarf being strongly affected by tidal inter- 
actions with the MW still now, although D M constitutes 
90% o f its total mass. In the same line, IMartin et al.1 

(2007) , following deep observations in the r band, infer 
that Willman 1 may probably be surrounded by tidal 
tails. The authors give as a plausible scenario that the 
dwarf could have been significantly tidally stripped but 
only in the past, when the object was more luminous and 
massive. At that age, Willman 1 could have lost most of 
its outskirts, only the innermost regions surviving intact. 
This picture would then allow the two contradictory ar- 
guments to coexist, since at least a correct modeling for 
the core of the dwarf may be possible assuming this re- 
gion to be in dynamical equilibrium at present. 

A well established model for the DM dis tribution, and 
the o nly one so far applied to Willman 1 (jStrigari et al.1 
12008ft. is the Navarro-Fr enk- White (NFW) density pro- 
file (jNavarro et al.lll997ft : 



PNFwO) 



r 

Ps{ r< 



1 + — 



(2) 



where p s and r s are a typical scale density and radius 
respectively. The astrophysical factor can be written as: 



J(*o) 



1 

47T 



dX[p 2 (r)*B^(6)}, (3) 



where "Jo denoted the direction of the target. The first 
integral is performed over the spatial extension of the 
source, the second is performed over the line of sight 
variable A. The density is convoluted with the Gaussian 
function B$ r (9) describing the telescope angular reso- 
lution where = \& — \&o is the angular distance with 
respect to the center of the object. We remark that the 
integration of Equation[3]involves foreground (MW halo) 
and extraga lactic background whose contr ibutions can be 
substantial (Elsa esser fc Mannheiml[2005l ). 

At a distance of 38 kpc, the scale radius corresponds 
to an extension of 0.54° in the sky, which is well in- 
side the MAGIC field of view (~ 3.5°), but is rather 
extended compared with the telescope angular resolu- 
tion of 0.1°. This evidence is mitigated by the fact 
that the main emission still comes from the very core 
of the source, due to the very steep NFW profile at 
the center and the square density dependence. For this 
reason, we perfor med an analysis adapted for slightly 
extended sources (|Sitarek fc Mirzovanl 12008) . To com- 
pute the astrophysical factor for Willman 1, we sub- 
stitute r s = 0.18 kpc and p s = 4 x 10 8 M /kpc 3 
taken from IStrigari et al.1 (|2008ft into Equation [2] and 



by computing Equation [3] we obtain J(^o) ~ 3.5 x 
10 17 GeV 2 /cm 5 . This value is possibly the largest over 
the rest of dSphs. 

3.2. Particle Physics Factor 

In many mSUGRA models, the lightest SUSY parti- 
cle is one of the four neutralinos (\i 4) which are linear 
combinations of the Bino B°, the Wino W° and the two 

neutral higgsinos H®, H 5. Only five free parameters 
fully characterize mSUGRA: the scalar mass too and the 
gaugino mass TO1/2 defined at the unification scale, the 
trilincar scalar coupling Aq and the ratio tan (3 of the 
Higgs vacuum expectation values. In addition, one needs 
to define the sign of the Higgs mass parameter sign(/i). 
The typical DM 7-ray annihilation spectrum is a contin- 
uum with a sharp cutoff at the DM candidate mass and 
possibly with a bump at energies larger than 0.6 m x in 
case the IB is present (see Figure [1]). 

The particle physics factor can be written as a prod- 
uct of two terms, the first depending only on the DM 
candidate mass and cross section, and a second term, 
depending on the 7-ray spectra of annihilation, which 
must be integrated above the energy threshold Eq of the 
telescope: 



2 m 2 



S(E) dE , 



(4) 



where 



J xx, 



is the total averaged thermal cross section 



times the relative velocity of particles, m x is the DM 
particle mass, and the factor 2 takes into account that 
the neutralino annihilates with itself. The 7-ray annihi- 
lation spectrum is composed of different contributions: 
S{E) = J2i dN^/dE where dN*/dE is the spectrum of 
the ith annihilation mode. 

The mSUGRA parameter space is conventionally de- 
scribed in a mo © m 1 / 2 plane, after having fixed the 
other free parameters. Usually, four zones are identi- 
fied: the bulk region, with low too and TO1/2 and neu- 
tralino masses at around 100 GeV, the focus point where 
mo and the neutralino are more massive, the funnel re- 
gion where both mo and TO1/2 take large values, and the 
co- annihilation tail characterized by large TO1/2. A neu- 
tralino shows different annihilation modes depending on 
the location in this p lane. A representative set of bench- 
marks was defined by (Battagfi^E^D (|2TOl|,[200J). Here- 
after, we use a subset o f four slightly mod i fied Ba ttaglia 
models, as defi ned by IBringmann et ail (|2008bft using 
DarkSUSY 4.01 (|Gondolo et al.ll2004ft . which include the 
contribution of IB in the computation of the cross sec- 
tions and spectra: models I', J', K\ F* for the bulk, 
co-annihilation, funnel and focus point regions respec- 
tively. All the defining parameters, as well as the derived 
$ pp (> E ), are listed in Tableffl 

The total estimated flux due to DM annihilation com- 
puted according to Equation [T] is given in Table [5J These 
estimations must be taken with carefulness since there 
are some sources of uncertainty that may largely af- 
fect the values for the predicted flux: (1) our lack of 
knowledge of the DM density profile may change the 
astrophysical factor by more than 1 order of magni- 
tude; (2) the presence of DM substructures in the dwarf 
may enhance the 7-ray flux at least by a factor of 
2 — 3 according to N— body simulations (|Diemand et al.1 
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TABLE 1 

Definition of benchmark models as in Bringma nn et alI (l2008Bf) and computation of the particle physics factor: m 1( / 2 and 
mo [GeV] are the gaugino and scalar mass respectively defined at the unification SCALE; tan p is the ratio of the Higgs 

EXPECTATION VALUES; Aq [GeV] THE TRILINEAR COUPLING CONSTANT AND sign(p) THE SIGN OF THE HlGGS MASS; m x [GeV] IS THE 
NEUTRALINO MASS; (<TV XX ) [CM 3 S" 1 ] IS THE CROSS SECTION TIMES THE RELATIVE VELOCITY OF DM PARTICLES AND & PP (> 100) [CM 3 

GeV -2 S" 1 ] IS THE PARTICLE PHYSICS FACTOR ABOVE 100 GeV. 
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Fig. 1. — Diffe rential particle physic s factor for the benchmarks 
models as in IBringmaniieTa l. (2008b). Line gamma emissions are 
not included, since their contribution to the flux is almost negligi- 
bly 

I2007allbl: iKuhlen et al.ll2 008) or even up to 20 according 



Martinez et al. (2009). Substructures are in fact 



to 

pected to be present not only in Willman 1 but also in 
any other DM halo, since cold dark matter (CDM) ha- 
los are approximately self-similar until a c utoff scale mass 
which lies in the range of 10~ 4 -10~ 12 M Q (|Profumo et al.l 
2006); (3) the exclusion of the baryons in the modeliza- 
tion of the total density profile. However, the effect of the 
adiab atic compression (jPrada et all 12004 iGnedin et all 
|2004| ). although important for larger DM halos, will prob- 
ably play a marginal role in the case of Willman 1 , given 
its relatively low amount of baryons even in the central 
regions, where the effect is expected to be more impor- 
tant; (4) the in clusion of the Sommerfe ld effect, recently 
investigated bv lLattanzi fc Silk (2008). This effect pre- 
dicts a general enhancement S of the flux, in case of very 
low DM particle velocities v, proportional to the inverse 
of the particle velocity: S oc l/v. In addition, for par- 
ticular values of the DM mass and very small velocities, 
due to the presence of bound states, the Sommerfeld ef- 
fect would give rise to resonances, which enhance the 
flux of a factor S oc l/v 2 . Finally, we underline that our 
choice of benchmarks does not scan the complete param- 
eter space and different ncutralinos could have a larger 
expected flux. A deeper study of the parameter space is 
therefore very desirable. 

4. MAGIC DATA 

The MAGIC (Major Atmospheric Gamma Imaging 
Cherenkov) telescope is located on the Canary Island 
La Palma (2200 m asl, 28.45°N, 17.54°W). MAGIC is 



currently the largest IACT, having a 17 m diameter 
tessellated reflector dish. The faint Cherenkov light 
flashes produced by atmospheric showers initiated by 
VHE 7-rays in the top atmosphere are reflected into 
an ellipsoidal image in the focal plane of the telescope, 
where a camera consisting of 577 photomultipliers (pix- 
els) records the ima ge. More details can be found in 
ICortina et all (|2005h . 

MAGIC observed Willman 1 between March and May 
2008. The source was surveyed at zenith angles between 
22° and 30°, which guarantees the lowest energy thresh- 
old. The source was tracked for 16.8 hr plus another 
9.3 hr in OFF observation mode, i.e. pointing to a dark 
patch in the sky close to Willman 1 where no 7-ray emis- 
sion is expected, for background estimation. The main 
background Cherenkov telescopes have to deal with is 
produced by cosmic hadronic particles impinging onto 
the top atmosphere and generating electromagnetic sub- 
showers that can mimic pure 7-rays showers, and by the 
night sky background. Background events are partly re- 
jected at the trigger level and in the off-line analysis event 
selection, following a procedure called gamma/hadron 
(g/h) separation. 

The analysi s proceeds as follow s (for a detailed de- 
scription, see lAlbert et al.l l2008bf ) : data are calibrated 
and the number of photoelectrons per pixel extracted 
([Albert et al.ll200"8cl ) , then an image cleaning selects pix- 
els with at least six photoelectrons (three photoelectrons 
in the boundary of the image) . Additional suppression of 
pixels containing noise is achieved by requesting a narrow 
time coincidence between adjacent pixels (~ 7 ns). Base d 
on the Hillas parameterization algorithm (Hillas 1985), 
the shower parameters are reconstructed. The hadronic 
background is supp ressed with a multivariate method , 
the Random Forest (|Breimanll200lHAlbert et al.ll2008dD . 
that uses the Hillas parameters to define an estimator 
called hadronness by comparison with Monte Carlo (MC) 
7-ray simulations. The hadronness expresses the likeness 
of an event to be a hadron and runs from for gammas to 
1 for hadrons. The Random Forest method is also used 
to estimate the energy of a reconstructed shower and the 
energy threshold is defined by the peak of the distribution 
of reconstructed MC gamma events. The g/h separation 
is optimized on a real data sample from the Crab Neb- 
ula, a supernova remnant and one of the brightest and 
stable 7-ray emitters, which is taken as standard candle 
in very high energy 7-ray astronomy. The optimization 
yields a best set of cuts in the Hillas parameters which 
defines the gamma and hadron acceptance of the analy- 
sis. In our case, the optimal set of cuts is obtained for 
an energy threshold of 100 GeV and a hadronness cut of 
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0.15. The overall data quality is very high, with only 7% 
data rejection, resulting in 15.5 hours effective observa- 
tion time. Independent cross-checks were performed on 
the data giving compatible results. 

5. RESULTS AND DISCUSSION 

No significant 7-ray excess beyond 100 GeV above the 
background was observed in 15.5 hr of observation of 
the sky region around Willman 1, according to Figure [21 
where the "a-plot" is reported. The a-parameter is the 
angular distance between the shower image main axis and 
the line connecting the image barycenter and the camera 
center. Due to their isotropic origin, hadronic events, 
in case they survive the analysis, are randomly oriented 
in the camera both in the ON-data and OFF-data sam- 
ple. This is reflected into a rather smooth distribution of 
events in the a-plot, the nonperfect flatness being due to 
an increased camera acceptance for showers with small 
a. On the other hand, 7-rays trace back the source, and 
thus the orientation of the shower image is toward the 
center of the camera. Therefore, in case of positive de- 
tection, an excess of events in the ON-data above the 
OFF-data sample is expected at small a. A fiducial re- 
gion a < 12° is chosen where the signal is assumed with 
a cut slightly larger than for a pointlike source to take 
into account the moderate source extension. The OFF- 
data are normalized to the ON-data in the region where 
no signal is expected, i.e., between a = 30° and a = 80°. 
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Fig. 2. — Willman 1 a-plot as seen by MAGIC in 15.5 hours above 
a fiducial energy threshold of 100 GeV and using a hadronness< 
0.15. The red crosses represent the ON-data sample, the blue 
shaded region is the OFF-data sample normalized to the ON-data 
sample between 30° — 80° . The vertical red dotted line represents 
the fiducial region a < 12° where the signal is expected. 

The si g nifica nce is calculated using Equation 17 of 
iLi fc Mai (|1983h . The number of excess events N exc {> 
100 GeV) = -223 ± 223 is calculated as the difference 
between the number of ON-events and the number of 
OFF-events i n the fiducial a— region. The method from 
iRolke etldl ((2005) was applied to estimate the upper 
limit in the number of excess events with a 90% confi- 
dence level and including 30% of systematic errors, giving 
as a result N£ c (> 100 GeV) - 191.4 events. This value 
is used to reconstruct the corresponding photon flux for 
a general 7-ray spectrum S(E) using: 



>£<o 



N, 



U . I . 



JS(E)A°ff(E) dEAT J Eo 



S(E)dE , (5) 



where Eq GeV is the energy threshold, Afff^E) is the ef- 
fective telescope area and AT is the effective observation 
time. 

We applied Equation [5] for the four neutralino bench- 
marks defined. Results are reported in Tabled where we 
also compare < J > "'^ Q with the estimated flux of Table [TJ 
We also report the upper limits on the boost factors that 
are required to match the two fluxes, again calculated 
separately for each neutralino model. The boost factor 
is defined as the ratio between the upper limit and the 
theoretical flux, and defines the minimal boost that the 
theoretical flux should be subject to in order to allow 
for a positive detection of the source. In order to pro- 
vide results less dependent of the particular benchmark 
spectrum, we also calculated flux upper limits in four 
different energy bins [100-170, 170-350, 350-1000, 1000- 
20 000] GeV for a generic annihilation spectrum without 
cutoff and spectral index —1.5. Respectively, the result- 
ing upper limits are [9.94, 4.75, 0.68, 0.35 ]xl0~ 12 pho- 
tons cm -2 s _1 . 

Table [2] reveals that although we derived upper lim- 
its of the same order of magnitude for the four models 
considered, there are evident differences in the prospects 
of detection for different neutralinos. The boost factor 
largely depends on the benchmarks, but the main differ- 
ences are connected to the particle physics factor, which 
varies by orders of magnitude among the different bench- 
marks, as shown in Table [T] The best prospects are for 
neutralinos in the funnel region (model K') of the pa- 
rameter space, for which the mass is large enough to place 
the cutoff well within the MAGIC energy threshold but 
still small enough not to reduce the particle physics fac- 
tor $ pp of Equation Q] too much. Next, with similar 
boost requirements follow the I' and F* models. In the 
former case, the effect of the IB plays an important role 
at energies close to the cutoff even if the neutralino mass 
is very close to the MAGIC energy threshold. In the lat- 
ter case, although the IB effect is negligible, the signal is 
very extended in the energy region suitable for MAGIC, 
whereas the large mass makes the flux suppression too 
large. The worst case scenario is the co-annihilation neu- 
tralino. In this case, even if the IB contribution is large, 
the intrinsic total cross section is very low, which makes 
the flux very low compared with the others. The IB effect 
cannot counteract this intrinsic deficit. 

Although the results of Table [2] are at first glance not 
looking overly promising for indirect detection of the dis- 
cussed models in the immediate future, one should keep 
in mind that the estimation of the flux is model depen- 
dent and there are intrinsic effects that could boost up 
the flux for up to 1 order of magnitude if the effect of sub- 
structures, astrophysical uncertainties, and/or the inclu- 
sion of baryons is considered (as already discussed in Sec- 
tion [22]) an d for several orders of magnitude, if the DM 
particle velocity is small enough to give rise to the Som- 
merfeld effect. This could change com pletely the outlook, 
as also discussed in lPieri et all ([2009D . 

6. CONCLUSION 

In the context of DM searches, we have observed the 
Willman 1 dwarf galaxy with the MAGIC telescope for 
a total of 15.5 hr between March and May 2008. Will- 
man 1 represents one of the best DM-dominated systems 
known in the universe to search for DM at present, ac- 
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TABLE 2 

Comparison of estimated integral flux above 100 GeV using Equation [1] for the benchmarks models defined in Table [T] and 

THE UPPER LIMIT IN THE INTEGRAL FLUX $ uL ABOVE 100 GeV COMING FROM MAGIC DATA IN UNITS OF PHOTONS CM -2 S _1 . On THE 
RIGHTMOST COLUMN, THE CORRESPONDING UPPER LIMIT ON THE BOOST FACTOR B u L REQUIRED TO MATCH THE TWO FLUXES IS 

CALCULATED. 



BM 


$ molfe, (> 100 GeV) 


100 GeV ) 


B u.l. 


/' 


2.64x10""* 


9.87 x lO" 1 ^ 


3.7 x 10 4 


J' 


4.29xl0~ 17 


5.69 x 10~ 12 


1.3 x 10 5 


K' 


2.32xl0~ 15 


6.83 x 10~ 12 


2.9 x 10 3 


F* 


2.09X10" 16 


7.13 x 10" 12 


3.4 x 10 4 



cording to its inferred dynamical properties and distance. 
No 7-ray signal was detected above an energy thresh- 
old of 100 GeV. We have obtained different flux upper 
limits of the order of 10~ 12 ph cm~ 2 s _1 separately for 
four benchmark models considered in the framework of 
mSUGRA. 

Using the latest estimations of its structural param- 
eters to build the DM density profile as well as the in- 
clusion of the recently proposed IB mechanism, we calcu- 
lated the boost factors needed to match the expected flux 
values from DM annihilation in Willman 1 with the upper 
limits obtained from the data. We can see that boosts 
in flux in the order of 10 3 are required in the most opti- 
mistic scenario considered. However, uncertainties in the 
DM distribution, the role of DM substructure or/and the 
contribution from the Sommerfeld enhancement may re- 
duce this required boost significantly. It is expected that 



deeper observations of the Willman 1 dSph with the up- 
coming MAGIC II telescope will allow us to improve the 
flux limits presented here by a factor 2 — 10. 
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